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Experiments demonstrate that the null space of the rigidity matrix determines grain 

motion during vibration-induced compaction 
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Using a previously developed experimental method to reduce friction in mechanically stable pack¬ 
ings of disks, we find that frictional packings form tree-like structures of geometrical families that lie 
on reduced dimensional manifolds in configuration space. Each branch of the tree begins at a point 
in configuration space with an isostatic number of contacts and spreads out to sequentially higher 
dimensional manifolds as the number of contacts are reduced. We find that gravitational deposition 
of disks produces an initially under-coordinated packing stabilized by friction on a high-dimensional 
manifold. Using short vibration bursts to reduce friction, we compact the system through many 
stable configurations with increasing contact number and decreasing dimensionality until the system 
reaches an isostatic frictionless state. We find that this progression can be understood as the system 
moving through the null-space of the rigidity matrix defined by the interparticle contact network in 
the direction of the gravitational force. We suggest that this formalism can also be used to explain 
the evolution of frictional packings under other forcing conditions. 

PACS numbers: 45.70.-n,61.43.-j,64.70.ps,83.80.Fg 


One of the great successes of statistical mechanics is 
the ability to calculate the average physical properties of 
macroscopic systems in thermal equilibrium, e.g. atomic 
and molecular liquids. The microstates and probabili¬ 
ties with which they occur are required for determining 
the average properties. For example, in the microcanon- 
ical ensemble, each of the microstates at fixed energy is 
equally probable. 

Here, we are interested in determining the structural 
and mechanical properties of static packings of frictional 
granular materials, which are collections of macroscale 
grains that interact via repulsive contact forces [l|. Our 
studies focus on frictional packings, which have not re¬ 
ceived as much attention as idealized packings of fric¬ 
tionless particles [2[ . Granular packings are also strongly 
out-of-equilibrium since thermal fluctuations do not give 
rise to particle rearrangements. Grain motion is induced 
instead through applied external loads Q. Moreover, the 
current state of the system may depend strongly on the 
protocol used to generate it In these systems, the 

definition of the appropriate microstates, their probabili¬ 
ties, and the framework for calculating average quantities 
is still under development M- 

In previous numerical simulations [l2| and experimen- 
tal studies 0, we showed that static packings of fric¬ 
tionless disks occur as distinct points in configuration 
space, whose number grows with the number of particles 
N in the system. We characterized each static packing 
using an invariant of the interparticle contact matrix and 
showed that static packings were not sampled with equal 
probabilities when generated using typical isotropic com¬ 
pression algorithms [T2| . 


In recent simulations, we showed that static packings 
of frictional disks do not occur as points in configuration 
space, but instead as lines, areas, and higher dimensional 
structures (z.e., geometrical families [l^) with a dimen¬ 
sion m that grows with the difference in the number of 
contacts of the system from the isostatic value for 
frictionless disks, m = — Nc > 0 [TqI. Thus, at 

each m, there are an infinite number of packings of fric- 
tionai disks, but they can be classified into a finite set of 
geometrical families with the same contact network. 

Granular systems can undergo rearrangements from 
one static packing to another in the presence of an ap¬ 
plied load, e.g., slow compaction during vertical tap¬ 
ping 0. There have been a number of theoretical 
studies of slow compaction in vibrated granular packings 
including Tetris-like and other lattice mode ls fl^ and 
adsorption-desorption or parking-lot models [l^, which 
have provided phenomenological understanding of the 
logarithmic relaxation of the density. Here, we seek a 
particle-level understanding of rearrangements that oc¬ 
cur during vibration in terms of the motion along ge¬ 
ometrical families and from one geometrical family to 
another. 

In this Letter, we describe quasi-two dimensional ex¬ 
periments of vibrated stainless steel particles and detailed 
analyses of the vibration-induced particle motion. We 
successively applied small amplitude, high frequency vi¬ 
bration bursts to the packings to break interparticle con¬ 
tacts and induce small-scale particle rearrangements. Us¬ 
ing accurate particle tracking techniques, we monitored 
changes in the interparticle contact network as the sys¬ 
tem evolved from dilute to more dense packings during 


2 




FIG. 1: (Color online) Sketch of the quasi-two dimensional 
vibration cell that contains a packing of seven bidisperse thin 
cylinders (4 small and 3 large with diameter ratio 1.246) with 
11 contacts. The amplitude of the vertical vibrations is in¬ 
dicated by A and gravity is directed downward. Rem and fn 
locate the center of mass of the packing and position of the 
nth particle. 


the vibrations. We found two key results: 1) The pack¬ 
ings that occur during the compaction process are not 
randomly distributed in configuration space, instead they 
form connected lines or ‘geometrical farnily trees’, which 
verifies our previous simulation results [l^. 2) The mo¬ 
tion of the grains during the compaction process occurs 
in the null space of an under-coordinated rigid link net¬ 
work El formed by the interparticle contacts. In par¬ 
ticular, the particle motion can be accurately predicted 
by projecting the direction of the gravitational force onto 
the null space. 

We performed quasi-two dimensional experiments of 
vibrated granular materials as shown in Fig. [H The vi¬ 
bration cell has width = 5.33 cm, height > 5 cm, 
and thickness Lz = 0.34 cm and contained N = 7 bidis¬ 
perse stainless steel cylinders (4 small and 3 large with 
diameters as = 1.26 cm and ai = 1.57 cm, masses 
rris = 2.8 g and mi = 4.3 g, and thickness 0.32 cm). 
We focused on systems with a small number of grains so 
that we can enumerate most of the geometrical family 
trees during compaction. The vibration cell is formed by 
two plastic side walls, two glass plates, and a movable 
piston on the bottom. The diameter of each cylinder 
is a factor of 4 larger than its thickness, which ensures 
that the cylinder axes are always parallel to each other 
and perpendicular to the direction of gravity. The piston 
is connected to a computer controlled electromagnetic 
shaker that vibrates the bottom boundary. 

For the experimental protocol, we first randomize the 
initial positions of the cylinders by applying high ampli¬ 
tude and low frequency (50 Hz) vibrations for one second 
and allow them to settle to a frictionally stabilized pack¬ 
ing. We then periodically apply small amplitude, high 


frequency (440 Hz) vibration bursts, each with duration 
10 ms. The acceleration from these bursts causes the 
contacts between the cylinders (and between the cylin¬ 
ders and walls) to break, which temporarily removes the 
frictional forces between particles. After each burst, the 
system settles under gravity, contacts reform, and the 
system forms a new static packing. Over a series of bursts 
the system compacts and the number of contacts grows 
on average, terminating on a ‘frictionless’ packing with 
Nc = = 2N = 14 contacts for systems with fixed 

boundaries and gravity. We analyzed 6901 initial pack¬ 
ings without rattler particles on the bottom boundary 
and applied 20 vibration bursts to each. After roughly 
15 bursts, these systems converged to 979 distinct fric¬ 
tionless packings. 

We imaged the system with a digital camera using 
backlighting. We recorded images of the system af¬ 
ter waiting 0.5 s following each burst. We tracked the 
positions of the cylinders with a spatial resolution of 
6 X 10“^(Js in both the horizontal and vertical direc¬ 
tions [l8|. 

In Fig. [2] (a), we show the evolution of the center 
of mass Rem and contact number for a set of packings 
that compacts to a particular frictionless packing with 
contacts. Mechanically stable frictional packings 
with fixed boundaries and gravity can occur in the range 
(3A^ - l)/2 < Nc < 2N or 10 < A^c < 14 for 

N = 7. The system is initialized in a dilute packing with 
Nc = 11 < N^^^. The center of mass evolves smoothly 
until the packing becomes mechanically unstable, and a 
new contact network is formed, in this case with Nc = 12. 
The Nc = 12 contact network smoothly evolves until it 
becomes unstable, and a new interparticle contact forms 
with Nc = 13. The evolution of the center of mass along 
a geometrical family is continuous as long as the con¬ 
tact network remains the same. Changes in the contact 
network are signaled by abrupt changes in the direction 
of motion of the center of mass. Since we apply small 
amplitude vibrations, the evolution of the center of mass 
terminates on a frictionless packing with contacts. 

When the system is initialized in a different packing, it 
can evolve to the same frictionless packing in Fig. [2] (a) 
or another one. In Fig. [2] (b), we show the geometrical 
family tree that leads to the frictionless packing shown in 
Fig. [2] (a). Note that the center of mass approaches that 
of the frictionless packing only along particular directions 
in configuration space. In Fig.[2](c), we show the centers 
of mass of all packings found in the experiments. Note 
that the most probable packings tend to be the most 
compact. 

We now describe a method to quantitatively predict 
the particle motion that occurs during the compaction 
process. The center of mass motion is characterized 
by smooth evolution while the contact network is fixed, 
punctuated by abrupt changes in direction when the con¬ 
tact network changes. While the contact network is fixed, 
we model the system by a network of rigid links between 
contacting particles. Since — Nc = m > 0^ the 
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FIG. 2: (Color online) (a) Evolution of the center of mass 
(Xcm,Ecm) as the system compacts and the contact number 
increases from Nc = 11 to 14. The star indicates a frictionless 
packing with = 14 contacts. Snapshots of this ‘geomet¬ 
rical family tree’ are included for each value of Nc. The solid 
line shows the solution of Eq. [J] in the overdamped limit for 
each of the contact networks, (b) Evolution of the centers of 
mass for all family trees (different line types) that terminate 
on the frictionless packing in (a), (c) All family trees found in 
the vibration experiments. The family trees with probability 
greater than 0.6% are shaded black, while the others appear 
gray. The dashed rectangular region indicates the area shown 
in (b). 


2N X 2N rigidity matrix K will possess m floppy eigen- 
modes e* = ..., ^In} form the null 

space and 2N — m non-floppy modes [ 13 . Each contact 
network will possess a different rigidity matrix and set of 
null-space modes. We define the rigidity matrix and its 
eigenmodes e* ^ for each experimental packing, where the 


subscripts s and n indicate the initial condition and vi¬ 
bration burst number, respectively. We have shown pre¬ 
viously that vibrations release the frictional contacts [Ij . 
We now hypothesize that the release of friction induces 
motion in the null space along directions that do not cost 
energy, and thus frictional packings live in the null space. 

To test whether the particle motion is confined to the 
null space, we calculate the normalized difference be¬ 
tween the total particle motion and the motion that oc¬ 
curs in the null space. 


W. 


(Ar^,n)^ - (Ar^,n ■ 

(Ar5,„)2 


( 1 ) 


where Ars,„ = rs,„+i - and is the 2N- 
dimensional vector that gives the positions of the par¬ 
ticles for initial condition s and burst n. Ws,n = 0 in¬ 
dicates that all particle motion is confined to the null 
space. In Fig. [3] (a), we plot the cumulative distribution 
C{Ws,n)^ which shows that 99% of the particle motions 
satisfy Ws,n ^ 10“^. Thus, nearly all of the particle 
motion during compaction is confined to the null space. 

To complete our theoretical analysis of vibration- 
induced compaction, we further hypothesize that particle 
motion occurs in the direction of gravity projected onto 
the null space. We now compare the direction of motion 
in configuration space during compaction to the direc¬ 
tion of the gravitational force decomposed onto to the 
null space. 


m 


i=l 


( 2 ) 


where g = {0, ..., 0, is a 2Wdimensional vec¬ 

tor with a zero for the x-component and the gravitational 
force for the ^-component for each particle. The degree 
to which the particle displacements track the direction 
of gravity can be determined by calculating the normal¬ 
ized overlap between Gs,n and To measure the 

deviation from complete overlap, we define 


= 1 - 


• Ar.s 


Gs^n 11 Ar^^T^ 


(3) 


as^n = 0 indicates that a given displacement Ar^^^ is 
aligned with the direction of gravity in the null space. In 
Fig. [3] (b), we plot the cumulative distribution C{as^n)^ 
which shows that 99% of the particle motions satisfy 
(^s,n ^ 10“^. Thus, nearly all of the particle displace¬ 
ments are in the direction of gravity in the null space. 
However, a few displacements out of 10^ possess as^n ^ 1 
(with overlap angles greater than 30°) indicated by the 
small peak in the probability distribution P{as,n) in the 
inset to Fig. [3] (b). These few instances of large as^n 
are caused by two effects: 1) The vibration bursts were 
not large enough to break some of the frictional contacts 
between particles and the bottom wall. In which case, 
these particles did not move in the direction of gravity 
in the null space. These occurrences of large as^n can 
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FIG. 3: (Color online) (a) Cumulative distribution C(Ws,n) 
of the normalized differences Ws,n between the particle mo¬ 
tion and the motion that occurs in the null space for each 
initial condition s and vibration burst n. (See Eq. [Tl) The 
open red circle indicates that 99% of the projections satisfy 
Ws,n ^ 10“^. The inset shows the probability distribution 
T^(Ws,n)- The minimum Ws,n ~ 10“^ is set by the resolution 
of the particle tracking, (b) Cumulative distribution C{as,n) 
of the deviations from complete overlap as,n between the par¬ 
ticle displacements and direction of gravity for each s and n. 
The open red circle indicates that 99% of the overlaps satisfy 
OLs,n ^ 10“^. The inset shows the probability distribution 
P{as,n). 

be removed by increasing the amplitude of the vibration 
or by adding constraints to the rigidity matrix to rep¬ 
resent frictional contacts. 2) The particle positions in 
the packing are such that different subsequent particle 
motions are equally likely. For example, when a parti¬ 
cle falls vertically downward on top of another particle. 
Whether the particle moves downward and to the right 
or to the left is sensitive to the precise horizontal loca¬ 
tion of the top particle. Once the direction of motion in 
configuration space is chosen, ag^n is nearly zero during 
subsequent evolution. 

We also solved for the full particle trajectories using 


the instantaneous rigidity matrices obtained from exper¬ 
iments by numerically solving Newton’s equations of mo¬ 
tion for the configuration space vector r: 

(Pv (jr 

where M is the mass matrix and b is the damping pa¬ 
rameter. We solved Eq. |4]in the overdamped limit, z.e., 
dr/dt = —Krib — gjb^ and find that the simulations are 
able to accurately recapitulate the full particle trajecto¬ 
ries in the experiments as shown by the solid line in Fig. [2] 
(a). 

We performed quasi-two dimensional experiments of 
granular materials undergoing vibration-induced com¬ 
paction. Using accurate particle tracking techniques, we 
monitored the motion of all grains in the system fol¬ 
lowing successive vibration bursts. Using the center of 
mass to classify the static packings, we find that they 
are organized into geometrical families. Abrupt changes 
in the direction of the evolution of the family in con¬ 
figuration space signal changes in the contact network. 
These experimental results represent a breakthrough in 
our understanding frictional granular packings and con¬ 
firm recent simulations that identified geometrical fam¬ 
ilies in isotropically compressed frictional packings [ 13 . 
We also modeled each contact network using rigid links 
between contacting grains and showed that grain mo¬ 
tion between vibration bursts occurs in the null space 
of the under-coordinated rigidity matrix. More specifi¬ 
cally, we find that the particles move in the direction of 
gravity projected into the null space. This novel method 
for determining particle motion can be applied to other 
driving mechanisms, e.g. fixed applied shear stress, and 
will allow us to understand shear-induced jamming and 
protocol-dependent mechanical properties in frictional 
packings. 
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